What, How, When, and Where: Spatiotemporal Water Quality Hazards of Cyanotoxins in Subtropical Eutrophic Reservoirs.

Though toxins produced during harmful blooms of cyanobacteria present diverse risks to public health and the environment, surface water quality surveillance of cyanobacterial toxins is inconsistent, spatiotemporally limited, and routinely relies on ELISA kits to estimate total microcystins (MCs) in surface waters. Here, we employed liquid chromatography tandem mass spectrometry to examine common cyanotoxins, including five microcystins, three anatoxins, nodularin, cylindrospermopsin, and saxitoxin in 20 subtropical reservoirs spatially distributed across a pronounced annual rainfall gradient. Probabilistic environmental hazard analyses identified whether water quality values for cyanotoxins were exceeded and if these exceedances varied spatiotemporally. MC-LR was the most common congener detected, but it was not consistently observed with other toxins, including MC-YR, which was detected at the highest concentrations during spring with many observations above the California human recreation guideline (800 ng/L). Cylindrospermopsin was also quantitated in 40% of eutrophic reservoirs; these detections did not exceed a US Environmental Protection Agency swimming/advisory level (15,000 ng/L). Our observations have implications for routine water quality monitoring practices, which traditionally use ELISA kits to estimate MC levels and often limit collection of surface samples during summer months near reservoir impoundments, and further indicate that spatiotemporal surveillance efforts are necessary to understand cyanotoxins risks when harmful cyanobacteria blooms occur throughout the year.

Concentration of total microcystins associates with nitrate and nitrite, and may disrupt the nitrogen cycle, in warm-monomictic lakes of the southcentral United States.

Cyanobacterial blooms and the toxins they produce pose a growing threat worldwide. Mitigation of such events has primarily focused on phosphorus management and has largely neglected the role of nitrogen. Previous bloom research and proposed management strategies have primarily focused on temperate, dimictic lakes, and less on warm-monomictic systems like those at subtropical latitudes. The in-lake conditions, concentration of total microcystins, and microbial functioning of twenty warm-monomictic lakes in the southcentral United States were explored in the spring and summer of 2021. Our data revealed widespread microcystins in lakes across this region, some of which exceeded regulatory limits. Microcystins were higher in the spring compared to the summer, indicating that warm-monomictic lakes, even across a large range of precipitation, do not follow the trends of temperate dimictic lakes. Microcystins were found in surface waters and bottom waters well below the photic zone, reflecting the persistence of these toxins in the environment. Principal components analyses showed a strong association between microcystins, nitrate + nitrite, and Planktothrix relative abundance and transcriptional activity. Many systems exhibited stronger denitrification in the spring, perhaps contributing to the decreased toxin concentrations in the summer. Counter to most sampled lakes, one lake with the highest concentration of total microcystins indicated nitrogen cycle disruption, including inhibited denitrification. These findings are relevant to mitigating cyanobacterial blooms and toxin production in warm-monomictic systems, and suggests a need to consider nitrogen, and not solely phosphorus, in nutrient management discussions.

Using YY supermales to destabilize invasive fish populations.

A plausible biocontrol strategy for the eradication of invasive species involves augmenting wild populations with genetically modified supermales. Supermales contain double YY chromosomes. When they are augmented into a wild population, destabilization and eventual extinction occurs over time due to a strongly skewed gender ratio towards males. Here, we employ a mathematical model that considers an Allee effect, but we have discovered through simulation that the presence of supermales leads to an increase in the minimal number of females needed for survival at a value higher than the mathematically defined Allee effect. Using this effect, we focus our research on exploring the sensitivity of the optimized supply rate of supermale fish to the initial gender ratio and density of the wild populations. We find that the eradication strategy with optimized supply rate of supermales can be determined with knowledge of reproductive rate and survival fitness of supermale fish.

Everything is not everywhere: can marine compartments shape phytoplankton assemblages?

The idea that 'everything is everywhere, but the environment selects' has been seminal in microbial biogeography, and marine phytoplankton is one of the prototypical groups used to illustrate this. The typical argument has been that phytoplankton is ubiquitous, but that distinct assemblages form under environmental selection. It is well established that phytoplankton assemblages vary considerably between coastal ecosystems. However, the relative roles of compartmentalization of regional seas and site-specific environmental conditions in shaping assemblage structures have not been specifically examined. We collected data from coastal embayments that fall within two different water compartments within the same regional sea and are characterized by highly localized environmental pressures. We used principal coordinates of neighbour matrices (PCNM) and asymmetric eigenvector maps (AEM) models to partition the effects that spatial structures, environmental conditions and their overlap had on the variation in assemblage composition. Our models explained a high percentage of variation in assemblage composition (59-65%) and showed that spatial structure consistent with marine compartmentalization played a more important role than local environmental conditions. At least during the study period, surface currents connecting sites within the two compartments failed to generate sufficient dispersal to offset the impact of differences due to compartmentalization. In other words, our findings suggest that, even for a prototypical cosmopolitan group, everything is not everywhere.

Resisting annihilation: relationships between functional trait dissimilarity, assemblage competitive power and allelopathy.

Allelopathic species can alter biodiversity. Using simulated assemblages that are characterised by neutrality, lumpy coexistence and intransitivity, we explore relationships between within-assemblage competitive dissimilarities and resistance to allelopathic species. An emergent behaviour from our models is that assemblages are more resistant to allelopathy when members strongly compete exploitatively (high competitive power). We found that neutral assemblages were the most vulnerable to allelopathic species, followed by lumpy and then by intransitive assemblages. We find support for our modeling in real-world time-series data from eight lakes of varied morphometry and trophic state. Our analysis of this data shows that a lake's history of allelopathic phytoplankton species biovolume density and dominance is related to the number of species clusters occurring in the plankton assemblages of those lakes, an emergent trend similar to that of our modeling. We suggest that an assemblage's competitive power determines its allelopathy resistance.

Interplay between r- and K-strategists leads to phytoplankton underyielding under pulsed resource supply.

Fluctuations in nutrient ratios over seasonal scales in aquatic ecosystems can result in overyielding, a condition arising when complementary life-history traits of coexisting phytoplankton species enables more complete use of resources. However, when nutrient concentrations fluctuate under short-period pulsed resource supply, the role of complementarity is less understood. We explore this using the framework of Resource Saturation Limitation Theory (r-strategists vs. K-strategists) to interpret findings from laboratory experiments. For these experiments, we isolated dominant species from a natural assemblage, stabilized to a state of coexistence in the laboratory and determined life-history traits for each species, important to categorize its competition strategy. Then, using monocultures we determined maximum biomass density under pulsed resource supply. These same conditions of resource supply were used with polycultures comprised of combinations of the isolated species. Our focal species were consistent of either r- or K-strategies and the biomass production achieved in monocultures depended on their efficiency to convert resources to biomass. For these species, the K-strategists were less efficient resource users. This affected biomass production in polycultures, which were characteristic of underyielding. In polycultures, K-strategists sequestered more resources than the r-strategists. This likely occurred because the intermittent periods of nutrient limitation that would have occurred just prior to the next nutrient supply pulse would have favored the K-strategists, leading to overall less efficient use of resources by the polyculture. This study provides evidence that fluctuation in resource concentrations resulting from pulsed resource supplies in aquatic ecosystems can result in phytoplankton assemblages' underyielding.

Lumpy species coexistence arises robustly in fluctuating resource environments.

The effect of life-history traits on resource competition outcomes is well understood in the context of a constant resource supply. However, almost all natural systems are subject to fluctuations of resources driven by cyclical processes such as seasonality and tidal hydrology. To understand community composition, it is therefore imperative to study the impact of resource fluctuations on interspecies competition. We adapted a well-established resource-competition model to show that fluctuations in inflow concentrations of two limiting resources lead to the survival of species in clumps along the trait axis, consistent with observations of "lumpy coexistence" [Scheffer M, van Nes EH (2006) Proc Natl Acad Sci USA 103:6230-6235]. A complex dynamic pattern in the available ambient resources arose very early in the self-organization process and dictated the locations of clumps along the trait axis by creating niches that promoted the growth of species with specific traits. This dynamic pattern emerged as the combined result of fluctuations in the inflow of resources and their consumption by the most competitive species that accumulated the bulk of biomass early in assemblage organization. Clumps emerged robustly across a range of periodicities, phase differences, and amplitudes. Given the ubiquity in the real world of asynchronous fluctuations of limiting resources, our findings imply that assemblage organization in clumps should be a common feature in nature.

Towards an Understanding of the Interactions between Freshwater Inflows and Phytoplankton Communities in a Subtropical Estuary in the Gulf of Mexico.

Subtropical estuaries worldwide face increased pressure on their ecosystem health and services due to increasing human population growth and associated land use/land cover changes, expansion of ports, and climate change. We investigated freshwater inflows (river discharge) and the physico-chemical characteristics of Galveston Bay (Texas, USA) as mechanisms driving variability in phytoplankton biomass and community composition between February 2008 and December 2009. Results of multivariate analyses (hierarchical cluster analysis, PERMANOVA, Mantel test, and nMDS ordination coupled to environmental vector fitting) revealed that temporal and spatial differences in phytoplankton community structure correlate to differences in hydrographic and water quality parameters. Spatially, phytoplankton biomass and community composition responded to nutrient loading from the San Jacinto River in the northwest region of the bay (consistent with nutrient limitation) while hydraulic displacement (and perhaps other processes) resulted in overall lower biomass in the Trinity River delta (northeast region). The influence of inflows on phytoplankton diminished along a north to south gradient in the bay. Temporally, temperature and variables associated with freshwater inflow (discharge volume, salinity, inorganic nitrogen and phosphorus concentrations) were major influences on phytoplankton dynamics. Dissolved inorganic nitrogen: phosphorus (DIN:DIP) ratios suggest that phytoplankton communities will be predominately nitrogen limited. Diatoms dominated during periods of moderate to high freshwater inflows in winter/spring and were more abundant in the upper bay while cyanobacteria dominated during summer/fall when inflow was low. Given the differential influences of freshwater inflow on the phytoplankton communities of Galveston Bay, alterations upstream (magnitude, timing, frequency) will likely have a profound effect on downstream ecological processes and corresponding ecosystem services.

Phytoplankton succession in recurrently fluctuating environments.

Coastal marine systems are affected by seasonal variations in biogeochemical and physical processes, sometimes leading to alternating periods of reproductive growth limitation within an annual cycle. Transitions between these periods can be sudden or gradual. Human activities, such as reservoir construction and interbasin water transfers, influence these processes and can affect the type of transition between resource loading conditions. How such human activities might influence phytoplankton succession is largely unknown. Here, we employ a multispecies, multi-nutrient model to explore how nutrient loading switching mode might affect phytoplankton succession. The model is based on the Monod-relationship, predicting an instantaneous reproductive growth rate from ambient inorganic nutrient concentrations whereas the limiting nutrient at any given time was determined by Liebig's Law of the Minimum. When these relationships are combined with population loss factors, such as hydraulic displacement of cells associated with inflows, a characterization of a species' niche can be achieved through application of the R* conceptual model, thus enabling an ecological interpretation of modeling results. We found that the mode of reversal in resource supply concentrations had a profound effect. When resource supply reversals were sudden, as expected in systems influenced by pulsed inflows or wind-driven mixing events, phytoplankton were characterized by alternating succession dynamics, a phenomenon documented in inland water bodies of temperate latitudes. When resource supply reversals were gradual, as expected in systems influenced by seasonally developing wet and dry seasons, or annually occurring periods of upwelling, phytoplankton dynamics were characterized by mirror-image succession patterns. This phenomenon has not been reported previously in plankton systems but has been observed in some terrestrial plant systems. These findings suggest that a transition from alternating to "mirror-image" succession patterns might arise with continued coastal zone development, with crucial implications for ecosystems dependent on time-sensitive processes, e.g., spawning events and migration patterns.

Phytoplankton assemblage characteristics in recurrently fluctuating environments.

Annual variations in biogeochemical and physical processes can lead to nutrient variability and seasonal patterns in phytoplankton productivity and assemblage structure. In many coastal systems river inflow and water exchange with the ocean varies seasonally, and alternating periods can arise where the nutrient most limiting to phytoplankton growth switches. Transitions between these alternating periods can be sudden or gradual and this depends on human activities, such as reservoir construction and interbasin water transfers. How such activities might influence phytoplankton assemblages is largely unknown. Here, we employed a multispecies, multi-nutrient model to explore how nutrient loading switching mode might affect characteristics of phytoplankton assemblages. The model is based on the Monod-relationship, predicting an instantaneous growth rate from ambient inorganic nutrient concentrations whereas the limiting nutrient at any given time was determined by Liebig's Law of the Minimum. Our simulated phytoplankton assemblages self-organized from species rich pools over a 15-year period, and only the surviving species were considered as assemblage members. Using the model, we explored the interactive effects of complementarity level in trait trade-offs within phytoplankton assemblages and the amount of noise in the resource supply concentrations. We found that the effect of shift from a sudden resource supply transition to a gradual one, as observed in systems impacted by watershed development, was dependent on the level of complementarity. In the extremes, phytoplankton species richness and relative overyielding increased when complementarity was lowest, and phytoplankton biomass increased greatly when complementarity was highest. For low-complementarity simulations, the persistence of poorer-performing phytoplankton species of intermediate R*s led to higher richness and relative overyielding. For high-complementarity simulations, the formation of phytoplankton species clusters and niche compression enabled higher biomass accumulation. Our findings suggest that an understanding of factors influencing the emergence of life history traits important to complementarity is necessary to predict the impact of watershed development on phytoplankton productivity and assemblage structure.

Ammonium treatments to suppress toxic blooms of Prymnesium parvum in a subtropical lake of semi-arid climate: results from in situ mesocosm experiments.

Prymnesium parvum is a haptophyte alga that forms toxic, fish-killing blooms in a variety of brackish coastal and inland waters. Its abundance and toxicity are suppressed by ammonium additions in laboratory cultures and aquaculture ponds. In a cove of a large reservoir (Lake Granbury, Texas, USA) with recurring, seasonal blooms of P. parvum, ammonium additions were tested in mesocosm enclosures for their ability to suppress blooms and their effects on non-target planktonic organisms. One experiment occurred prior to the peak abundance of a P. parvum bloom in the cove, and one encompassed the peak abundance and decline of the bloom. During 21-day experiments, weekly doses raised ammonium concentrations by either 10 or 40 µM. The added ammonium accumulated in experimental mesocosms, with little uptake by biota or other losses. Effects of ammonium additions generally increased over the course of the experiments. The higher ammonium dose suppressed the abundance and toxicity of P. parvum. The biomass of non-haptophyte algae was stimulated by ammonium additions, while positive, negative and neutral effects on zooplankton taxa were observed. Low ammonium additions insufficient to control P. parvum exacerbated its harmful effects. Our results indicate a potential for mitigating blooms of P. parvum with sufficient additions of ammonium to coves of larger lakes. However, factors excluded from mesocosms, such as dilution of ammonium by water exchange and sediment ammonium uptake, could reduce the effectiveness of such additions, and they would entail a risk of eutrophication from the added nitrogen.

A chemical approach for the mitigation of Prymnesium parvum blooms.

Known as Golden Algae in popular media, the harmful algal bloom causing organism Prymnesium parvum secretes increased amounts of toxic chemicals called prymnesins when stressed, resulting in major fish kills in Texas. Although many options exist for mitigation of blooms, a feasible protocol for control of blooms on large-scale impoundments has yet to be identified. Chemical control of P. parvum using six different enzyme inhibiting aquatic herbicides was explored in laboratory experiments. Of the six chemicals screened, one (flumioxazin) was selected for further study due to a significant decrease in P. parvum cell numbers with increasing chemical concentration. It was applied to natural plankton communities during in-situ experiments (Lake Granbury, Texas). The first experiment was conducted during a period of P. parvum bloom initiation (March) and the second experiment conducted during a post bloom period (April). Experiments were carried out in 20 L polycarbonate carboys covered in 30% shade cloth to simulate natural light, temperature and turbulence conditions. Through cell counts via light-microscopy, the chemical flumioxazin was found to cause significant decreases in P. parvum, but no significant differences in zooplankton abundance during the period of bloom initiation. However, significant decreases in adult copepods were observed during the post bloom period, with no significant decreases in P. parvum most likely due to decreased light penetration and inhibition of the photosensitive mode of action.

STIMULATING EFFECT OF ANABAENA SP. (CYANOBACTERIA) EXUDATE ON PRYMNESIUM PARVUM (HAPTOPHYTA)(1).

Prymnesium parvum blooms have become more frequent in the south-central United States, leading to significant ecological and economic impacts. Allelopathic effects from cyanobacteria were suggested as a mechanism that might limit the development of P. parvum blooms. This research focused on the effects of cultured cyanobacteria, Anabaena sp., on P. parvum. Over a 6-d period, daily additions of filtrate from the senescent Anabaena culture were made to P. parvum cultures growing in log phase. All treatments, including several types of controls, showed reductions in P. parvum biomass over the course of the experiment, but the treatments receiving Anabaena filtrate were reduced to a lesser degree, suggesting that filtrate from the senescent cyanobacteria culture was beneficial to P. parvum in some way. This unexpected outcome may have resulted from stimulation of heterotrophic bacteria by the addition of Anabaena filtrate, which likely contained exudates rich in dissolved organic carbon compounds. P. parvum was then able to supplement its nutritional requirements for growth by feeding on the elevated bacteria population. These findings coupled to previous observations suggest that interactions between cyanobacteria and P. parvum in natural environments are complex, where both allelopathic and growth-stimulating interactions are possible.

Enrichment and detection of Escherichia coli O157:H7 from water samples using an antibody modified microfluidic chip.

Low abundant (<100 cells mL(-1)) E. coli O157:H7 cells were isolated and enriched from environmental water samples using a microfluidic chip. The poly(methylmethacrylate), PMMA, chip contained 8 devices, each equipped with 16 curvilinear high aspect ratio channels that were covalently decorated with polyclonal anti-O157 antibodies (pAb) and could search for rare cells through a pAb mediated process. The chip could process independently 8 different samples or one sample using 8 different parallel inputs to increase volume processing throughput. After cell enrichment, cells were released and enumerated using benchtop real-time quantitative polymerase chain reaction (PCR), targeting genes which effectively discriminated the O157:H7 serotype from other nonpathogenic bacteria. The recovery of target cells from water samples was determined to be approximately 72%, and the limit-of-detection was found to be 6 colony forming units (cfu) using the slt1 gene as a reporter. We subsequently performed analysis of lake and wastewater samples. The simplicity in manufacturing and ease of operation makes this device attractive for the selection of pathogenic species from a variety of water supplies suspected of containing bacterial pathogens at extremely low frequencies.

A mechanistic explanation for pH-dependent ambient aquatic toxicity of Prymnesium parvum carter.

The harmful algal bloom species Prymnesium parvum has caused millions of dollars in damage to fisheries around the world. These fish kills have been attributed to P. parvum releasing a mixture of toxins in the water. The characterized toxins, reported as prymnesin-1 and -2, have structural similarities consistent with other known ionizable compounds (e.g., ammonia). We investigated whether pH affects the toxicity of P. parvum under conditions representative of inland Texas reservoirs experiencing ambient toxicity from bloom formation. We evaluated pH influences on toxicity in laboratory and field samples, and modeled the physicochemical properties of prymnesins. Aquatic toxicity to a model fish and cladoceran was reduced by lowering pH in samples obtained from reservoirs experiencing P. parvum blooms; similar observations were confirmed for experiments with laboratory cultures. A pKa value of 8.9 was predicted for the prymnesins, which suggests that ionization states of these toxins may change appreciably over surface water pH of inland waters. These findings indicate that ionization states of toxins released by P. parvum may strongly influence site-specific toxicity and subsequent impacts to fisheries. Consequently, these results emphasize the importance of understanding processes that affect pH during P. parvum blooms, which may improve predictions of ambient toxicity.

Mixing of supersaturated assemblages and the precipitous loss of species.

Mechanisms influencing species richness are many. Recent theoretical research revealed additional mechanisms that involved neutral and lumpy coexistence and alternating assemblage states. These mechanisms can lead to conditions where the number of coexisting species is greater than the number of limiting resources, that is, species supersaturation. Our research focused on the role of disturbances (migration and pulsed through-flows) in supersaturated plankton systems. Our simulations employed 30 different supersaturated assemblages generated by using various ecological principals. Our findings indicated that immigration rates as low as 0.1% of total biomass per day generally led to regional homogenization of species and dramatic extinction events, with assemblages characteristic of lumpy coexistence being more resilient than those characteristic of neutral coexistence or alternating states. Generally, pulsed through-flows tended to offset, to some extent, the negative effects of migration. The precipitous loss of species with the onset of migration is observed in other systems as well, for example, cichlid fish communities of East Africa rift lakes and songbird assemblages from Indian Ocean islands. While many explanations have been offered to explain postimmigration extinctions in species-rich systems, another explanation might be that the assemblages in these systems are in a fragile state of supersaturated coexistence.

Laboratory tests of ammonium and barley straw extract as agents to suppress abundance of the harmful alga Prymnesium parvum and its toxicity to fish.

Prymnesium parvum is a harmful alga whose blooms can cause fish kills in brackish waters. Two potential suppressants of this alga were tested, ammonium and barley straw extract (BSE), at temperatures of 10, 20 and 30 degrees C. Laboratory batch cultures were grown for 3 weeks at each temperature, with weekly doses of ammonium or BSE at either low or high levels, or a no-dose control treatment. The growth rate of P. parvum during exponential phase was highest at 20 degrees C and lowest at 10 degrees C, and was stimulated by the highest ammonium dose. Only cultures grown at 20 degrees C were toxic to fish. The highest ammonium dose abolished such toxicity and reduced the endpoint population density of P. parvum. BSE did not reduce the exponential growth rate, endpoint density, or toxicity to fish of P. parvum. The results support the use of ammonium additions, but not BSE, to suppress harmful blooms of P. parvum in those circumstances where the possible disadvantages can be managed.